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1. Introduction 

Nuclear and relativistic effects in atomic spectra are treated in the pioneer works of 
Stone [I], [2j and Veseth ||. The theory of the mass shift has then been reformulated 
by Palmer 0. Calculations of nuclear motional effects in many-electron atoms have 
been performed by Parpia and co-workers f| |6| in the relativistic scheme, using fully 
relativistic wave functions, but adopting the non-relativistic form of the recoil operator. 
Relativistic nuclear recoil corrections to the energy levels of multicharged ions have been 
estimated by Shabaev and Artemyev |7j who derived the relativistic corrections of the 
recoil Hamiltonian. In a study of isotope shifts of forbidden transitions in Be- and B-like 
argon ions, Tupitsyn et al |§ showed that a proper evaluation of the mass isotope shift 
requires the use of this relativistic recoil operator. The latter has also been shown to be 
crucial by Porsev et al |J for calculating isotope shifts of transitions between the fine 
structure energy levels of the ground multiplets of Fe I and Fe II. 

As far as computational atomic structure is concerned, the extension of the available 



relativistic codes such as grasp2k |TI| or MCDF-gme [11], is needed for estimating 
these mass corrections properly for any many-electron system. Programs to calculate 
pure angular momentum coefficients for any scalar one- and two- particle operator are 



available [ 1 3 1 but do require the knowledge of the tensorial structure of the operators 



to be integrated between the many-electron atomic wave functions [113. The tensorial 



form of the nuclear recoil Hamiltonian is derived in the present work, opening interesting 
perspectives for calculating isotope shifts in the multiconfiguration Dirac-Hartree-Fock 
(MCDHF) framework. 

2. The relativistic mass shift operator 

In the MCDHF method, the atomic state function (ASF) ^(-jPJMj), of a stationary 
state of an atom, is expressed as a linear combination of symmetry-adapted configuration 
state functions (CSFs) $>(-y p PJMj), i.e. 

m^PJMj) = Y^cMlpPJMj), (1) 

p 

where J is the total electronic angular momentum of the state, 7 represents the electronic 
configuration and intermediate quantum numbers, and P stands for the parity. The 
mixing coefficients c p and the one-electron radial wave functions spanning the CSFs are 
optimized by solving the MCDHF equations iteratively until self-consistency. The latter 
are derived by applying the variational principle to the energy functional based on the 



Dirac-Coulomb Hamiltonian 14 



N N 1 

£ (c Qi • Pl + fa - l)c 2 + V( n )) +J2~ ( 2 ) 
. . ^ / r^j 



N N 

Hdc ^_ 

=1 i<j 
where V(rj) is the monopole part of the electron-nucleus interaction, a and /3 are the 
(4 x 4) Dirac matrices and c is the speed of light (c = 1/a in atomic units, where a is 
the fine-structure constant). 
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The mass shift of the energy levels in an atom with nuclear mass M is caused by 
the recoil motion of the atomic nucleus. The corresponding recoil Hamiltonian 

1 V"^ / OlZ ( (cKj • 7*j) Tj\ \ 

has been derived within the lowest-order relativistic approximation and to first order in 
m/M by Shabaev and collaborators ||. Rewriting it as the sum of the normal mass 
shift (NMS) and specific mass shift (SMS) contributions and using the tensorial form 
r 1 = rC 1 , (fD becomes 

Hms = H^ms + Hsms , (4) 



with 



#nms = ^ E (Pi ~ ■Vi-°^ 7 {pi i - C]) C] ■ , (5) 

#sms = [Pi-Pj-^ a i-Pj-^-( a i- C }) C l-Pj) , ( 6 ) 

i<j v 1 1 / 

that, in both cases, are rewritten as a sum of three separate contributions: 

#NMS = -f^NMS + -^NMS + -^NMS > (?) 

and 

Hsms = #sms + -^sms + h sms • ( 8 ) 
Since the expectation values of the NMS and SMS operators are evaluated with the 
MCDHF wave functions, the expectation values (H^ MS ) and (Hg us ) partly contain 
the relativistic contributions. Tupitsyn et al || pointed out that averaging the non- 
relativistic recoil operator with the relativistic wave functions strongly overestimates 
the relativistic correction to the recoil effect such as it becomes important to use the 
complete form (Q) when one works in the relativistic scheme. 

2.1. Normal mass shift expectation value 

The (mass- independent) normal mass shift parameter -Knms is defined by the following 
expression 

= (V(lPJMj)\HKus\y(lPJMj)) . (9) 



-K'nms 
M 



3 

NMS- 



By analogy with (0), we define i^NMS as the sum of -R'nms = -^nms + -^nms + K 
Applying the Wigner-Eckart theorem jn|, the matrix element of the normal mass shift 
operator is Mj-invariant and is proportional to the reduced matrix element (r.m.e.)H 

^ = ^= {mpjwvuswnipj)) 

= [*( 7 PJ)||# NMS ||*( 7 PJ)] . (10) 



I The two definitions of r.m.e are related to each other through (7J'||O fe |j7J) = \J2J' + 1 ['y«7 / ||0||'yj"]. 
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Using multiconfiguration expansions (|1|), the reduced matrix elements of the general 
spherical tensor operator T k becomes 

[v&( 7 PJ)||T fc ||vI/( 7 PJ)] = J2 [®(%PJ)\\T k \m ls PJ)] . (11) 

p,s 

The reduced matrix elements of the one-electron operator T k = '^2 i t k (i) between CSFs 
is expressed as a sum over single-particle reduced matrix elements 

m lp PJ)\\T k \m ls PJ)} = T ps( ab ) [n a K a \\t k \\n b K b ] , (12) 
where the T ps (ab) are the spin-angular coefficients arising from Racah's algebra [|l|, 



14] ) 16 . Introducing the one-body normal mass shift operator associated to 



(H~NMS — ^NMS(^) 

hnus = ± (p 2 ~ (« + (« • Cl ) Cl ) ■ p) . (13) 

we hereafter derive the expression of its r.m.e., using relativistic central- field one-electron 
wave functions 

r \tQn a , Ka {r)tt- Ka ,rn a {0,<j),O-) J 

P a and Q a are respectively the large and small components of the relativistic one-electron 
radial wave function a = (n a K a ), where k = (I — j)(2j + 1). 

Introducing the notation d r = 4-, the action of the operator p 2 on the large (F = P) 
and the small (F = Q) component of a relativistic wave function 

p 2 ^1q k ^9, h&) = \ (-d 2 r + F ntK (r)n K>m (9, 0, a) (15) 

is found using 



1 i 2 

p 2 = -A = -- 2 d T r 2 d r + - . (16) 

From this expression and integrating by parts, the first term of the one-electron matrix 
element NMS operator ( |i~3D becomes, 

p 2 

(n a K a m a \ — \n b K b m b ) = 5(K a m a , n b m b ) (17) 
1 [°° j f , a r,^, a , a „,, a „^ , h(h + l)P a Pb + h(! b + l)QaQb 



*\ j dr((d r P a )(d r P b ) + (d r Q a )(d r Q b ) + 
with I = 2j — I. Building the Dirac matrices a from ct = o~ x <g> a with 

one rewrites the second and the third parts of the NMS operator ( pT3|) as 

(^)(a.p+(a.C)(C.p)) S «r.®A=(° M (19) 
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-ctZ 
~ 7 h r 



(cT-p+itT-C 1 ) (C l -p)). 



Taking into account that (see (A. 4. 9) and (3.2.14) in fT4|| ) 

cr ■ C 1 — a • e r = oy ; C 1 ■ p = e r ■ p = (—id r ) 

and 

<T ■ p = —io r o r H 

\ r 

with K = — (1 + <t • I), the operator A becomes 

, -aZ\ , x / „ K + l 
A = I — ] (-ia r ) 2d r + 



2r J " ' \ r 
Acting on the one-electron relativistic wave function component, it gives 



A- 



F„.Jr) 



-a,.. 



-aZ\ i ( „ k — 1 . 

2d r + ] F n Jr)n„ K , 



(20) 

(21) 
(22) 

(23) 
(24) 



v \ 1r I r \ r 

from which one derives, integrating by parts, the one-electron matrix element of the 
second and third parts of the NMS operator 

(n a K a m a \a x <g> A\n b n b m b ) = 5(K a m a , K b m b ) (25) 

x ij~ ^ qjm+qap. + { _ aZ) ^ (w> + 

Combining (|T7|) and ( p^|) to deduce the r.m.e. of the normal mass shift operator (|T3|), 
we obtain the final expression 

[n « ||/i NMS ||n 6 « 6 ] = 5(n a ,K, b ) (26) 



x 



2M Jo 



(d r P a ) (d r P b ) + (d r Q a ) (d r Q b ) + 



l b (l b +l)P a P b + l b (l b + l)Q a Q b 



n ry\ Qad r Pb + Qbd r P a . , ry\ ( K b ~ ^- \ . , . n n , ) 

-2aZ) h (-qZ) — — [Q a Pb + W^) K"*. 



Specific mass shift expectation value 

Similarly to fllOD, the (mass- independent) specific mass shift parameter K$ms is defined 
as 

-Ksms 



M 



[*( 7 PJ)||^sms||*(7^] , 



(27) 



and K$ MS , -^sms' -^sms as ^ s contributions according to (§). Its evaluation requires the 
calculation of the corresponding matrix elements in the CSF space. For the general 
scalar two-particle operator 

G = Y,9<i,i) ( 28 ) 



with 



i<j 

9(hj) = Y,9^i,r j ) (T fc (.) • T k (j)) 
k 



(29) 
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the reduction of the many-electron r.m.e. in terms of the two-electron integrals X k , also 
called effective interaction strengths [ |T4|j , 

m lp PJ)\\J29(iJ)\MlsPJ)] = J2 Y. v{ vs^cd)X k {abcd) (30) 



i<j 



abed k 



can be performed using Racah's algebra [13|, [L4], [16| . For the specific mass shift 
Hamiltonian (^), using k — 1, all three terms have the particular form 

g{i,j) = gi{r i )g l {r j ){T\i).T 1 {3)) (31) 

in which the radial part gi(ri,rj) of ( f29"|) is factorized. Adopting the covariant notation 
for the 3j-symbol of Wigner and using the definition of the scalar product of two 
irreducible tensor operators and the Wigner-Eckart theorem, the matrix element of (j31|) 
can be written as follows 

i 

(ab\g(i,j)\cd) = 

q=-l 

where 



m a 1 j c 
3a q rn c 



m b q j d 
jb 1 rn d 



X\abcd), (32) 



X\abcd) = -(aWg^T'^ibWgir^md) . 
From the structure of @, the latter has three components 

X l (abcd) = Xl(abcd) + X^(abcd) + Xl(abcd) 
that we analyzed hereafter separately. 

2.2.1. First part: X\(abcd) 

Building X\(abcd) from (|33"D , we have 

g{ ri )T l (i)=p 1 (i) ; g(r J )T 1 (j)=p\j) . 
Introducing the one-electron reduced matrix element 

(allp^c) = -i^aWC 1 ^)^ {n a K a ,n c K c ), 

where V (nn, ti'k') is the Vinti radial integral 

d k{k + 1) - k'{k' + 1) 



(33) 
(34) 



'o 



+ 



QnK \J) 



dr 
' d 

dr 



2r 



(35) 
(36) 

P n , K ,(r)dr (37) 



-k(-k + 1) + k'(-k' + 1) 
2r 



Q n 'K'(r)dr 



the first contribution to the effective interaction strength writes as 

Xl(abcd) = (K a ||C 1 ||K c )(K fe ||C 1 ||/t (i )^(?2 a K a ,n c fi; c )^(n & fi; ;) ,n d K d ) , (3£ 
recovering the uncorrected relativistic expression used in || H|. 
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2.2.2. Second part: X 2 (abed) 

For the second term of the SMS operator, we identify from (|J), (§), (0) and fl34|) 

g (r l )T 1 (t) = ^(a x ®* 1 ( i )) ; g(r j)T x (j) = p x (j) . (39) 
Introducing the matrix element 
(a|r _1 (a x ® crj) |c) 

— (-Qa-P C (-/«a"lak g |/«c"lc) + QcPa (K, a m a | | ~ K c m c )) , (40) 

and using the r.m.e. fl36|), the corresponding contribution to the effective interaction 
strength is 

Xl(abcd) = -(K b \\C 1 \\K d )v(n b K b ,n d K d ) 

(^^— J (-^(-/Calk 1 !!^) + QcPaKlk'H " K C » • (41) 

££5. Third part: X\(abcd) 

Similarly, the two components of the third term of the SMS operator are 

g(n)T\i) = ^ (a x ® (a r (i)C l (i))) ; g^T^j) = p\j) . (42) 

Using the matrix element 

(air" 1 (a x g> (cr r Cj)) |c) 



/ — (Q a P c (-K a m a |Cg| - n c m c ) - Q c P a (K a m a \C\\K c m c )) , (43) 
Jo r 

and the r.m.e. (|36|), the third contribution to the effective interaction strength takes the 
form 

X\(abcd) = - (KbllC 1 ^^) (K a \\C 1 \\K c )<V (n b K b ,n d K d ) 

x f°° d r(—\(Q a P c -Q c P a ) , (44) 



r 

where we take advantage of (— « a ||C || — k c ) = (K a \\C k \\K c ) 



2.3. Useful one- electron reduced matrix elements 



Equations (f26|), (pq) , ([41]) and (|44]) are the final key expressions of the relativistic mass 
shift one-electron r.m.e. that involve the following three reduced angular one-electron 
matrix elements 

Mc 1 !^) = (-iy a+1/2 VU^m I _\ c /2 y(i a ,i c ,i), (45) 

(-K a \\^\\K c ) =5(l,l c )(-l) T ^ 2+ ^ 1 VW^ul (j 2 ^ I ,(46) 
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3c 1/2 1 



where n (l a , l c , 1) is defined by : 



-K (I I 1) = { 1 if ^ + 1 + ^even, 
otherwise. 



3. Applications 

We wrote a new program, hereafter referred as RMS2, for estimating the expectation 
values of the relativistic nuclear recoil operators using MCDHF wave functions 
calculated with the GRASP2K package fl(|. This code is based on the previous program 
SMS92 || in which 

• for the NMS, the one-electron radial integrals (expression (39) of the original 
paper 0) are replaced by the corresponding relativistic expression fl26|) , 

• for the SMS, the first contribution X\(ahcd) (expression (40) of the original 
paper |J) is corrected by adding the relativistic contributions ( f41~D and (|44]) . 

It is important to notice that the program SMS92 calculates the uncorrected 
NMS as the expectation value (%2iTi), where Tj is the Dirac kinetic energy operator 
Tj = c CKj ■ Pi + (A — 1 ) c 2 associated to electron i, while the program RMS2 
uses more accurately (H^ MS ) = (%2iPj/2M), which is consistent with section 2.1. An 
equivalent version has been written for the code MCDF-gme. 



In the present work, we evaluate the NMS and SMS parameters (§) and (|27 ) for 
some low-lying levels of neutral lithium, boron-like argon and two medium- carbon-like 
ions (Ca XV and Sc XVI) to investigate the importance of the relativistic corrections. 
The nuclear charge distribution is described by a Fermi model. Nuclear masses (Mjy) 
are calculated by taking away the mass of the electrons and the binding energy from 
the atomic mass (Ma), using the formula: 

M N (A 1 Z)=M A (A 1 Z)-Zm e + B a (Z) (49) 

where the total binding energy of the electrons (expressed in eV) is estimated 
using 



18, 19 



B el (Z) = 14.4381 Z 2m + 1.55468 • 1(T 6 Z 5 - 35 . (50) 
The atomic and nuclear masses relevant to the present work are reported in Tables 0. 
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Table 1. Atomic masses (Ma) [@ and nuclear masses (M^r) (in u) calculated from 
( ^9|) and ( |50| ) for lithium and argon isotopes. 



Isotope 


M A 


M N 


6 Li 


6.015122795(16) 


6.01386737 


7 Li 


7.01600455(8) 


7.01474907 


36 Ar 


35.967545106(29) 


35.9576862 


40 Ar 


39.9623831225(29) 


39.9525242 



When discussing a transition mass isotope shift, one needs to consider the variation 
of the mass parameter from one level to another. The line k frequency isotope shift, 
Su Al ' A2 = (SE Al ' A ' 2 — SE Al,A2 )/h, between the isotopes A 1 and A 2 , of nuclear masses 
Mi and M 2 respectively, is usually written as the sum of the normal mass shift (NMS), 
specific mass shift (SMS) and field shift (FS) contributions : 

OV k — Of fe NMS -h 01S k sus -\-0U k FS , (01) 



with 



AuM ( M 2 -M 1 \ AK MS ( M 2 -M x \ ~ 
^ = \ M\M~2~ ) ~h~ = \ M\M~2~ ) A ^ MS ' (52) 
where AKms is the difference of the Kms parameters of the levels involved in 
transition k. As far as conversion factors are concerned, we use^| Ai^Ms[GHz u] = 
3609.4824 Ai\" MS [m e .E h ] . Note that thanks to the separability enhanced in (f|), (|52]) can 
be applied to both the mass contributions NMS and SMS, separately. 



3.1. Hydrogen-like selenium 

Below we present some relevant calculations of the expression ( p6| ) for a heavy one- 
electron ion (Se XXXIV, Z = 34). This choice is motivated by the interesting 
comparison with the unpublished work of Kozlov ||22|| . The normal mass shift values 
calculated with the operators -^nms arLC1 (-^nms + -^nms)> usrn g the RMS2 program, are 
reported in Table ^|. In the second and third column respectively, comparison is made 
with the numerical results of Kozlov together with our analytical values. The latter are 
based on analytical hydrogenic wave functions^ [^|. The agreement is very satisfactory. 



§ This conversion factor is calculated as (m e /u)2R 00 c x 1.10" 9 = 3609.4824 using the 2006 CODATA 
recommended values of the fundamental physical constants |f2~lf . 

|| The values reported in Tables f| and || are based on or 1 = 137.035989500 adopted in GRASP2K. 
For lsi/2, the analytical result for -FCnms becomes 656.358886872 if adopting the or 1 — 137.035999679 
2006 CODATA value §3 • 
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Table 2. Contributions to the normal mass shift -KTnms parameters (in m e -Eh) for 
hydrogen-like selenium (Z = 34). 

RMS2 Kozlov 1 22j Analytic 



K 1 

-"-NMS 



ls 1/2 656.3589797 656.3589 656.358899684 

2p 1/2 154.8937900 154.8938 154.893789883 

2p 3/2 147.5192507 147.5192 147.519250700 

-^NMS + ^NMS 

lsi/2 -78.3588949 -78.3589 -78.3588996839 

2px/2 -8.0987869 -8.0988 -8.0987868710 

2p 3/2 -3.0192507 -3.0193 -3.0192506997 



3.2. Lithium-like systems using Dirac one-electron wave functions 

The SMS parameters for Li-like iron (Z = 26) and selenium [Z = 34) are calculated 
in the single configuration approximation using three-electron wave functions built on 
unscreened Dirac solutions. The results are reported in Table |^ and compared with 



independent estimations using an adapted version of MCDF-gme jlT| , |24] and with the 
analytical results. The three sets are consistent with each other but sensitively different 
from Kozlov's values |22| reported in the last column of the table. Note that the 



comparison is somewhat unfair to Kozlov since the grid parameters used for the discrete 
representation of orbital wave functions have been adapted in both programs (rms2 and 
MCDF-gme) to achieve a better accuracy. 
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Table 3. Contributions to the specific mass shift -K"sms ( m ™>eEh) parameters for 
Li-like iron (Z — 26) and selenium (Z — 34) using unscreened Dirac one-electron wave 
functions. 

RMS2 MCDF-gme Analytic Kozlov j22 



Li-like iron 






K 1 

^SMS 




ls 2 2 Pl/2 
ls 2 2p 3/2 


2 po 

r l/2 
2 po 

^3/2 


-55.24725067 
-53.26443137 


-55.24725061 
-53.26443136 


-55.247250683 
-53.264431362 


-55.2474 
-53.2645 


ls 2 2 Pl/2 
ls 2 2p 3/2 


2 po 

^1/2 
2 po 

*3/2 


3.48269304 
1.20278262 


-^SMS + 

3.48269307 
1.20278262 


^SMS 

3.482693070 
1.202782617 


3.4278 
1.1960 


Li-like selenium 




-"■SMS 




ls 2 2 Pl/2 
ls 2 2p 3/2 


2 po 

^1/2 
2 po 
^3/2 


-97.71464163 
-91.70637651 


-97.71464140 
-91.70637651 


-97.714641685 
-91.706376511 


-97.7150 
-91.7069 


ls 2 2 Pl/2 
ls 2 2p 3/2 


2 po 

^1/2 
2 po 

^3/2 


10.53884731 
3.55081371 


K 2 + 

SMS + 

10.5388474 
3.5508137 


^SMS 

10.53884746 
3.55081372 


10.2546 
3.5164 



3.3. Neutral lithium in the MCDHF approach 

The MCDHF active space method consists in writing the total wavefunction as a 
configuration state function expansion built on a set of active one-electron orbitals. To 
investigate the convergence of the property, the orbital set is systematically expanded 
up to n = 10, but imposing the angular restriction l max = 6 {i orbitals). The sequence 
of CSFs Active Spaces (AS) is resumed as follows 



AS 


= ls 2 2s, 


AS 2 


= AS + {2 P } , 


AS 3 


= AS 2 + {3s, 3p, 3d} , 


AS 4 


= AS 3 + {4s,4p,4d,4/}, 


AS 5 


= AS 4 + {5s, bp, bd, bf, bg} , 


AS 6 


= AS5 + {6s, 6p, 6d, 6f, 6g, 6h} , 


AS 7 


= AS 6 + {7s,7p, 7d, 7f, 7g,7h,7i}, 


AS 8 


= AS 7 + {8s, 8p, 8d, 8f, 8g, 8h, 81} , 


AS 9 


= AS 8 + {9s, 9p, 9d, 9f, 9g, 9h, 9i} , 


AS 10 


= AS 9 + {10s, lOp, lOd, 10 /, lOg, lOh, 10i} , 
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where the (nZ)-notation implies the relativistic shell structure j = I ± 1/2. The 
configuration space is increased progressively, by adding at each step a new layer of 
variational orbitals, keeping the previous ones frozen from the (n — 1) calculation. 
The MCDHF expansions are based on single and double (SD) excitations from the 
configuration reference. Triple excitations are investigated through SDT-configuration 
interaction (CI) calculations. 

Table 4. Uncorrected (if^ MS ) and corrected {K~nms) normal mass shift parameters 
(in m e Eh) for Li I. 



AS n SD SDT 







"NMS 


-P^NMS 


K 1 

NMS 


Knms 


n= 


■5 


ls 2 2s 2 , 
7.479955285 


5"l/2 

7.473188966 


7.480387512 


7.473620714 


n= 


■6 


7.480757179 


7.473989593 


7.481294538 


7.474526401 


n= 


■■7 


7.480843823 


7.474076167 


7.481413156 


7.474644913 


n= 


■■8 


7.482617678 


7.475849092 


7.483709525 


7.476940080 


n= 


9 


7.482764972 


7.475996085 


7.483865298 


7.477095534 


n= 


■10 


7.482767804 


7.475998981 


7.483872626 


7.477102933 


n= 


=5 


ls 2 2p 2 P° /2 
7.411878601 7.405201316 


7.412125687 


7.405448151 


n = 




7.412307495 


7.405629843 


7.412599631 


7.405921698 


n = 


=7 


7.412593434 


7.405916172 


7.413034981 


7.406357367 


n= 




7.414193990 


7.407516244 


7.415203718 


7.408525555 


n= 


9 


7.414351543 


7.407673608 


7.415377292 


7.408698936 


n= 


:10 


7.414365017 


7.407687009 


7.415402512 


7.408724081 


n = 


=5 


1, 2 2 P 2 P° /2 
7.411871260 7.405208436 


7.412118064 


7.405455006 


n= 




7.412300503 


7.405637317 


7.412592555 


7.405929108 


n= 


=7 


7.412584979 


7.405922413 


7.413026010 


7.406363146 


n= 


= 8 


7.414185599 


7.407522728 


7.415193271 


7.408530126 


n = 


= 9 


7.414343399 


7.407680348 


7.415366987 


7.408703663 


n = 


= 10 


7.414356793 


7.407693666 


7.415392260 


7.408728857 



Tables § and |5] present the evolution of the NMS and the SMS parameter, 
respectively. In each table both the uncorrected (K^ s ) and corrected (K MS ) values 
are reported. Comparing the SD and SDT calculations, we observe that the influence 
of the triple excitations reaches more than 1 % for the SMS while it is one order of 
magnitude smaller (0.1%) for the NMS. 

In Table |5| the individual contributions to the mass shift AK MS (= AK M g/h) 
parameters are reported for the 2p 1 / 2 2 P°/2 ~ 2s 2 Si/ 2 (Di line) and 2p 3 / 2 2 P^ 2 ~ 
2s 2 S\/2 (-D2 line) transitions in lithium. Values are calculated with the SD and SDT 
n = 10 active space final results of Tables [| and Although many robust theoretical 
studies on the resonance line transition isotope shifts are available (see Table [7| and 
discussion below), the comparison with other theoretical works presented in Table ^ is 
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Table 5. Uncorrected (i^g M g) and corrected (K$ms) specific mass shift parameters 
(in m e Eh) for Li I. 



AS n SD SDT 







K 1 

-"■SMS 




K 1 

-"-SMS 


-f^SMS 


n= 


5 


ls 2 2s 2 S 1/2 
0.3010343291 0.3008225633 


0.3013767853 


0.3011648528 


n= 


6 


0.3010361585 


0.3008243666 


0.3014579841 


0.3012459847 


n= 


7 


0.3019544951 


0.3017423943 


0.3024396569 


0.3022273237 




8 

O 


0.3018617523 


0.3016497153 


30241 1 5843 


3021 992791 


n= 


9 


0.3017987398 


0.3015867742 


0.3023512554 


0.3021390203 


n= 


10 


0.3018561821 


0.3016442119 


0.3024141615 


0.3022019200 


n= 


■5 


ls 2 2p 2 P° /2 
0.2489342617 0.2487564343 


0.2490604867 


0.2488826064 


n= 


6 


0.2482881614 


0.2481107089 


0.2484282378 


0.2482507247 


n= 


7 


0.2482993836 


0.2481222326 


0.2484107486 


0.2482335819 


n= 


Q 
O 


0.2476015693 


0.2474248242 


H OA 1 A KKTOOK 

U.Z4 ( 400 i ZZo 


H OA 70701 T A Q 


n= 


9 


0.2475207029 


0.2473440589 


0.2473719731 


0.2471954806 


n= 


10 


0.2476566450 


0.2474799659 


0.2475131224 


0.2473365865 


n= 


■5 


\s 2 2p 2 P° /2 
0.2489331884 0.2487377216 


0.2490592224 


0.2488636586 


n= 


6 


0.2482892038 


0.2480941432 


0.2484289962 


0.2482338261 


n= 


7 


0.2483022860 


0.2481070388 


0.2484139538 


0.2482185782 


n= 


8 


0.2476039924 


0.2474089620 


0.2474585594 


0.2472634414 


n= 


9 


0.2475233735 


0.2473284034 


0.2473750253 


0.2471799634 


n= 


10 


0.2476606363 


0.2474656360 


0.2475176569 


0.2473225644 



limited to the recent large-scale configuration-interaction Dirac-Fock-Sturm calculations 
of Kozhedub et al P3 since these authors precisely focus on the estimation of the 
relativistic nuclear recoil corrections. Kozhedub et al 's values are very consistent with 
our results: They report A(K^ MS + K^ MS ) = 0.33 and 0.38 GHz u for the D\ and D 2 
transitions respectively. However, the uncorrected NMS contribution and therefore, the 
total NMS values, sensitively differ from each other by around 1.6 GHz u. This latter 
discrepancy is not understood yet and clearly deserves further investigations. 

The uncorrected contribution of the SMS is also compared with the non-relativistic 
result of Godefroid et al |26| using the multiconfiguration Hartree-Fock method. More 
interesting is the comparison with the recent SMS values of Kozhedub et al [25| 
investigating the relativistic recoil corrections and using the same NMS and SMS 
partition according to (H) and (^). As for the NMS, the relativistic corrections are 
in very nice agreement (they report A(i^"g MS + -Kg MS ) = 0.12 and 0.06 GHz u for the 
Di and D 2 lines ) but the uncorrected forms do differ substantially with our estimation 
(they report for instance for the D x line, Ai^g MS = —198.78, against our value of 
-198.164 GHz u). 
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The comparison with observation for the individual mass contributions is also 
limited. There are a few reasons for this. First, as illustrated by (pi]), the field 
shift contribution should be properly subtracted from the observed transition frequency 
before trying to extract the mass contribution. But this is usually the other way 
round that makes the theoretical calculation of mass shifts interesting: for a few- 
electron atomic systems like lithium indeed, the difference between the mass contribution 
calculated by elaborate ab initio calculations and the observed transition IS allows to 
extract the change in the mean square charge radius of the nuclear charge distributions 
for all isotopes, as illustrated by the very recent and complete work of Nortershauser et 
al |27| . Another good reason is that once the FS "eliminated" , a clean separation of the 
NMS and SMS contributions could be criticized, as pointed out by Palmer However, 
remembering that for lithium, the FS is roughly 10 4 times smaller than the MS, it is 
worthwhile to neglect it for trying the mass separation exercise. There is indeed one 
experimental work by Radziemski et al p8| discussing the NMS and SMS separation 
in this line but as we will observe later (see Table 0), the corresponding experimental 
transition IS values are not aligned with most of the other observed values. In their 
work, these authors separate the two mass shift contributions from the experimental 
transition IS in 6 ' 7 Li, neglecting the field shift contribution and approximating the Bohr 
mass shift by the experimental observed level energy, as suggested by Martensson and 
Salomonson [p9|j , 

a bms = _^ e b -^E exp . (53) 
M M 

From the same expression, we build the transition Bohr mass shift for the 6 ' 7 Li isotope 
pair 

from the obserbed transition energy. Combining (|52|) and fl54|), one finds 



AKnms - rn e ^ E ^ xp = m e Uexp (55) 

from which we estimate the "observed" NMS values reported in Table || using 
the most recent absolute frequency measurements of Das and Natarajan []30[| . The 
corresponding "observed" A_K"sms values are calculated by substracting the so-estimated 
NMS contribution from the experimental IS line shifts (—443.9490(16) GHz u and 
—443.9126(29) GHz u, for D x and D 2 , respectively). Note that we did not take the 
liberty of reporting the frequency uncertainties estimated by Das and Natarajan on the 
separate contributions, the separability of NMS and SMS being by itself questionable. 

Cleaner and in principle less problematic should be the comparison of the total 
mass shifts, as reported in Table 0. On the theoretical side, we refer to the study 
of Korol and Kozlov |JT| treating electron correlation with configuration interaction 
(CI) and many-body perturbation theory (MBPT) methods with Dirac-Fock orbitals, 
to the calculations of Kozhedub et al f25f using large-scale configuration-interaction 



Dirac-Fock-Sturm method and to the Yan et al |32| calculations estimating the mass 
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Table 6. Individual contributions to the mass shift AKms (GHz u) parameters for 
the 2p 2 P°/ 2 ~ 2s 2 Si/2 and 2p 2 P%/ 2 ~ 2 5 1 / 2 transitions in lithium 









2po 2q , 
M/2 _ ^1/2 


2 po 2 o , 
*3/2 _ ^1/2 


NMS 


/KK 1 

NMS 


SD 


-246.899 


-246.928 






SDT 


-247.142 


-247.179 






SD 
SDT 


0.328 
0.333 


0.382 
0.384 




Ai^NMS 


SDT 


-246.812 


-246.795 




Other theory" 




-245.15 


-245.11 




Obs. fc 




-245.103 


-245.108 


SMS 


A-f^SMS 


SD 


-195.632 


-195.618 






SDT 


-198.164 


-198.148 




A(-^SMS + ^SMs) 


SD 
SDT 


0.127 
0.129 


0.061 
0.062 




Ai^SMS 


SDT 


-198.035 


-198.086 




Other theory" 




-198.78 


-198.77 




Other theory (NR) C 




-198.66 


-198.71 




Obs. d 




-198.843 


-198.101 



a CI Dirac-Fock-Sturm calculation of Kozhedub et al [|25|. 

b NMS values deduced from the transition frequencies [ |30| using ([35]) (see text). 
c Non- relativistic MCHF calculations p(| . 

d SMS values obtained by subtracting the "observed" NMS (see footnote b above) from 
the IS measured by Das and Natarajan raj. 



corrections from highly correlated non-relativistic wave functions expressed in Hylleraas 
coordinates^. From all these elaborate results, we only kept the mass contributions, 
systematically excluding the contributions from the nuclear size corrections. We already 
noticed the differences between Kozhedub et al 's results and ours appearing in the 
separate NMS and SMS contributions. As commented above, these differences do not 
arise from the relativistic corrections (K 2 + K 3 ), but rather from the "uncorrected" K 1 
values, and should be further investigated. Our results seem to be of higher quality than 

% The values of Yan et al reported in the Kozhedub's paper |25f suggest that the atomic mass has 
been used in order to evaluate the mass shift parameter. In table 0, the Yan et al 's values have been 
reevaluated using the nuclear mass. 
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the CI+MBPT results of Korol and Kozlov. As far as the differences with Yan et al 's 
results are concerned, we should keep in mind i) that our orbital active set is truncated 
to l max = 6, ii) that the layer approach adopted in the SD-MCDHF optimization could 
be a limiting factor and iii) that the convergence of the AKms parameter as a function 
of the size of the active set is slow and not yet achieved at n — 10, as illustrated by the 
comparison of the two n = 9 and n = 10 sets of results reported in the Table [7[ 

On the experimental side, we display in the same Table [7|, the experimental isotope 
shift values somewhat abusively converted in AT^ms parameters, ie. neglecting the FS 
contribution and inverting (|52|) , AK MS = 5i/ fc (M 1 M 2 )/(M 2 — Mi). As already mentioned, 
this conversion is unfair to physicists who do some huge efforts to extract the nuclear 
charge radii from the FS ||27|| , but has the merit of illustrating where the present modest 
contribution lies in the distribution of experimental values. From this not exhaustive 
chronological list ( |33| , j28| [34] , [35], |37|, [30]), it is clear that Radziemski et al 's results 
lie a bit outside the experimental distribution. 



Table 7. Mass shift AK M s (GHz u) for the 2p 1/2 2 P° /2 - 2s 2 S 1/2 and 2p 3/2 2 P° /2 - 
2s 2 Si/2 transitions in lithium, compared with experimental IS. 





2po 2 c , 
^1/2 _ ^1/2 


2po 2c , 
^3/2 ~ ^1/2 


Rcf. 


n = 9 


-445.1941 


-445.2330 


This work 


71 = 10 


-444.8442 


-444.8808 


This work 


Other theory 


-447(12) 


-447(12) 


Korol and Kozlov |3l| 




-443.81(20) 


-443.82(20) 


Kozhedub et al [g5| 




-443.860337(253) 


-443.876984(253) 


Yan et al || 


Experiment" 


-443.89(3) 


-443.91(2) 


Sansonetti et al |33|] 




-443.46(63) 


-443.59(63) 


Radziemski et al p8[ 




-443.9033(63) 


-443.9791(63) 


Scherf et al [B4[ 




-443.9045(29) 




Bushawet al [p5| 




-443.951(5) 




Walls et al [|36| 




-443.941(3) 


-443.948(4) 


Noble et al p7| 




-443.9490(16) 


-443.9126(29) 


Das and Natarajan |3(| 



a inverting (|||), ie. using AK M s = 6v(M 1 M 2 )/{M 2 - Mi) (see text) 



3.4- B-like argon 

Large-scale calculations are performed for ls 2 2s 2 2p 2 -P 1 °/ 2 3/2 °f B-like argon (Z = 18). 
The radial orbital basis is obtained from SD-MCDHF calculations, including single 
and double excitations from all shells of the {ls 2 2s 2 2p, ls 2 2p 3 } complex to increasing 
orbital active sets, up to the {10s9p8d7f6g3hli}. Subsequently to this layer-by- layer 
SD-MCDHF orbital optimization, RCI calculations are performed including the Breit 
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and QED effects in a space generated by SD excitations from the extended {ls 2 2s 2 2p, 
ls 2 2p 3 , Is 2 2s2p3d, ls 2 2p3d 2 } multireference set to the full orbital set. The expansion 
for the two J values includes more than 200 000 relativistic CSFs. This computational 
strategy has been developed by Rynkun et al []33 for the evaluation of transition rates 
in boron-like ions, from N III to Zn XXVI. 

Table ^ illustrates the convergence of the NMS and SMS contributions with the 
increasing of the active set. In Table §, the isotope shifts of the forbidden transitions 



Table 8. NMS and SMS parameters (in m e E h ) values for the states ls 2 2s 2 2p 2 P\ 
and ls 2 2s 2 2p 2 P° /2 states of B-like Ar. 



ASn ^NMS (-^NMS + ^NMs) ^SMS (^SMS + ^SMs) 



ls 2 2s 2 2p 2 P° 



n = 


3 


417.6959959 


-12.90905322 


-16.14960237 


0.46905612 


n = 


4 


418.1118745 


-12.91188239 


-16.30895958 


0.47763419 


n = 


5 


418.2177551 


-12.91178616 


-16.33759742 


0.47889811 


n = 


G 


418.2537558 


-12.91247514 


-16.32403619 


0.47580853 


n = 


7 


418.2919516 


-12.91360934 


-16.34490331 


0.47658869 


n = 


8 


418.2945374 


-12.91354179 


-16.34259843 


0.47612308 


n = 


9 


418.2974119 


-12.91373611 


-16.34355766 


0.47621031 


n = 


10 


418.2982181 


-12.91375324 


-16.34352453 


0.47609622 


n = 


10 ^expand 


418.2995162 


-12.91378955 


-16.33796596 


0.47600086 



ls 2 2s 2 2p 2 P° 



n = 


3 


417.3566053 


-12.66720414 


-15.91469860 


0.12821557 


n = 


4 


417.7946482 


-12.66774060 


-16.09016993 


0.13464620 


n = 


5 


417.9013634 


-12.66662806 


-16.11510903 


0.13335729 


n = 


6 


417.9378902 


-12.66759777 


-16.10142189 


0.13128456 


n = 


7 


417.9773669 


-12.66844094 


-16.12255333 


0.13135116 


n = 


8 


417.9799742 


-12.66842909 


-16.12002982 


0.13104864 


n = 


9 


417.9828424 


-12.66858843 


-16.12097800 


0.13104262 


n = 


10 


417.9836553 


-12.66862952 


-16.12094341 


0.13100231 


n = 


10 ^expand 


417.9846675 


-12.66865788 


-16.11555533 


0.13089443 



ls 2 2s 2 2p 2 P°/ 2 ~ 2 P%/2 m 36 ' 40 Ar are presented and compared with the results of 
Tupitsyn et al [§]. In their work, the CI Dirac-Fock method was used to solve the 
Dirac-Coulomb-Breit equation and to calculate the energies and the isotope shifts. The 
CSFs expansions were generated including "all single and double excitations and some 
part of triple excitations" . The nuclear charge distribution is described by a Fermi model 
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and is therefore consistent with the present work. The large cancellation of the terms 
involved in the transition isotope shift makes accurate calculations very challenging. 

Table |9] shows the individual contributions of operators (||) and (|J) to the 
wavenumber mass shift. A good agreement is observed between the two sets of values, 
the total wavenumber mass shift values differing by less than 0.8%. This example 
beautifully confirms the importance of the relativistic corrections to the recoil operator: 
the total wavenumber mass shift would be indeed 50% smaller if estimated from the 
uncorrected form of the mass Hamiltonian (H^ MS + Hg MS ). The nice agreement with 
Tupitsyn et al 's results is also a good sign of reliability for the tensorial form derivation 
of the nuclear recoil Hamiltonian of section 2. 



Table 9. Individual contributions to the wavenumber mass shift 8a (cm 1 ) for the 
forbidden transition ls 2 2s 2 2p 2 P°, 2 - 2 P° /2 in boron-like 3M0 Ar. 







Sa 












(^NMS + ^NMs) 






Total 


This work 


0.1054 


— 0.0821 


-0.0745 


0.1155 


0.0645 


Tupitsyn et al Q 


0.1053 


— 0.0822 


-0.0742 


0.1151 


0.0640 



3.5. C-like ions calculations 

As another illustration of the importance of the relativistic corrections to the recoil 
operator, the values of the SMS, NMS and total level mass shift parameters are 
reported in Table [H^ for the levels arising from the ground configuration ls 2 2s 2 2p 2 
in Ca XV and Sc XVI. As far as the calculations are concerned, the orbitals are 
obtained by SD-MCDHF calculations, considering single and double excitations from 
all shells of the {ls 2 2s 2 2p 2 , ls 2 2p 4 } Layzer complex to the {8s7p6d5 f Ag2h} active set. 
These MCDHF calculations are followed by relativistic configuration interaction (RCI) 
calculations, including the Breit interaction and the QED corrections, using the enlarged 
multireference {ls 2 2s 2 2p 2 , ls 2 2p 4 , ls 2 2s2p 2 3d, ls 2 2s 2 3<i 2 } set. The size of the expansions 
is around 350 000 relativistic CSFs. This computational method has been used by 



Jonsson et al [|39| to calculate transition rates, hyperfine structures and Lande g factors 
for all carbon-like ions between F IV and Ni XXIII. 

On the absolute scale of level shift parameters, one observes that the relativistic 
corrections (K^ MS + K^ MS ) to the NMS have the same order of magnitude than the 
uncorrected SMS contribution. Transition isotope shifts are more interesting properties 
since they are the real observables if the resolution is good enough. These are monitored 
by the differential effects on the level IS. It is interesting to infer from Table |1(] the 
possible mass isotope shifts on the intraconfiguration (M1/E2) transition frequencies. 
Considering for example the Ca XV 3 Pi — > 3 p2 transition, the uncorrected total mass 
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shift change is enlarged by a factor two when including the (K^s + ^ms) relativistic 
corrections. For the 3 Pi — > l D 2 transition, a similar increase of the mass shift is 
predicted but of "only" 20%. Some reduction could occur: this is the case of 3 P 2 — > 1 Sq 
(13%). For the 3 Pq — > 1 D 2 transition, the relativistic recoil corrections reach 48%. 



Table 10. Specific mass shift i^sMS, normal mass shift i^NMS , total mass shifts Kyis 
parameters (all in m e E\) for 2s 2 2p 2 levels of Ca XV and Sc XVI from multireference 
RCI calculations. K = K 1 + K 2 + K 3 . 

SMS NMS Total 
Level J Kg MS K SMS ^nms ^nms Kms k ms k ms + -^ms 

CaXV 



2s 2 2p 2 3 P -41.993 

1 -41.736 

2 -41.588 
2s 2 2p 2 1 D 2 -41.451 
2s 2 2p 2 1 S -41.862 



-40.659 558.358 

-40.741 558.070 

-40.821 557.854 

-40.756 557.479 

-41.199 557.000 



538.033 516.364 

537.983 516.334 

537.929 516.266 

537.607 516.028 

537.200 515.138 



497.374 -18.990 

497.242 -19.092 

497.108 -19.158 

496.851 -19.177 

496.001 -19.137 



Sc XVI 



3p 





-47.460 


-45.771 


621.221 


596.303 


573.761 


550.532 


-23.229 




1 


-47.117 


-45.874 


620.842 


596.236 


573.725 


550.362 


-23.363 




2 


-46.950 


-45.964 


620.592 


596.169 


573.642 


549.219 


-24.423 


l D 


2 


-46.765 


-45.922 


620.145 


595.826 


573.380 


549.904 


-23.476 







-47.229 


-46.421 


619.624 


595.390 


572.395 


548.969 


-23.426 



4. Conclusion and outlook 

The irreducible tensorial form of the nuclear recoil Hamiltonian is derived in the 
present work, opening interesting perspectives for calculating isotope shifts in the 
multiconfiguration Dirac-Hartree-Fock framework. We implemented the formalism in 
the relativistic package GRASP2K by writing a dedicated code (rms2) for estimating 
the expectation values of the relativistic nuclear recoil operators. The comparison 
with other works is satisfactory and the results are promising, although not achieving 
the accuracy of the state-of-the-art methodology available for a few-electron systems. 
Electron correlation remains the major problem that might be solved in our schema 
with the use of "localized pair-correlation functions interaction method", as proposed 



by Verdebout et al 1 40| . The present work enhances the fact that the relativistic 
corrections to the nuclear recoil are definitively necessary for getting reliable isotope 
shift calculations. The new computational tool, that we developed on the basis of the 
irreducible tensorial operator techniques, will hopefully provide valuable mass isotope 
shift data for large systems for which there are no reliable theoretical or experimental 
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